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ABSTRACT 

It is shown that there is a simpler way to derive the average guiding 
centre drift of a distribution of particles than via the so-called 
"single particle" analysis. Based on this derivation it is shown 
that the entire drift formalism can be considerably simplified, and 
that results for low order anisotropies are more generally valid than 
is usually appreciated. This drift analysis leads to a natural al= 
ternative derivation of the drift velocity along a neutral sheet. 

Full derivations are given in Burger, Moraal and Webb (1985; re= 
ferred to as BMW), of which copies will be available at the Conference. 


1. Alternative Derivation of Drift Velocity. A particle at position x, with 
momentum £, velocity v, in a magnetic field B has a gyroradius R = b x £/q, 
where b = B/B 2 . Its guiding centre is at Xg = A - £. Differentiation of 

x g w.r.t. time in the steady state (3/3t=0), with employment of the Lorentz 
force A = q(E. + v. x J3) then gives 


< v 


< 1 || > + Ixb + <£X (v.-V)b >/q, 


(D 


for the instantaneous, average drift velocity of particles inside a volume ele= 
ment d 3 x around x. For a particle distribution function F(x., £, t) with den= 
sity N, this average is defined as 


<• • •> 


= (1/N) Fd 3 p = (1/N) | 


dpp 2 


Fdfi, 


(2a) 


where dfi= sinOd0dc|). We shall also use differential (omnidirectional) 
averages in momentum space in the interval (p, p + dp) denoted by 




FdQ/ Fdn 


(2b) 


In terms of directions in momentum space, £ may be written as 


£ = p(cos0e x + s i n 0 cos <|> e. 2 + sinOsm<t>e. 3 ), (3) 

with 0 the pitch and 4> the phase angle relative to B. In index notation the 
last term of (1) is e ijk <p j v^>b k £ /q and determination of < v^ > or < v g > u 

only requires finding the components < p^ > of the < £V> tensor, given a par= 

ticular F. In BMW it is shown that the most complicated F which gives a 

tractable expression for v is 

9 co 

F(x,£,t)=Fgg+{F 1 QCOS0+F 11 sinOcos<()+F 1 ^_ 1 sin0sin(j)}+ £ F^gP^(cos0), (4) 


where P^ are the Legendre polinomials. 


<J) v > = 


<p. v > 

K ii ii 


0 0 
i<P_ L V ± > 0 
0 i<p v > 


The < £ y. > tensor is of the form: 

< p „ v „ > = i <pv( 1 + lm )> 

iK W >= < ^2 > = < *>3 W 3 > 


= 7 < pv (1 


1 F20\ 


(5) 
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For this form of F, the average drift velocity (1) becomes 

<V = <V + ixb4n* < P 1 V lxb + (<p l , v „ > ' (6) 


There are six noteworthy points about this seemingly familiar result: (i ) It 

is the instantaneous average guiding centre velocity of particles in the volume 
element d 3 x in terms of the fields £ and £ ^ns^de d 3 x. The usual "single par= 
tide" derivations (e.g. Rossi and Olbert, 1970) give a similar result, but for 
the time (gyroperiod) averaged drift velocity in terms of the fields Eg and E5g 
at the gutdtng centres. (n) Our analysis does not need the usual Single par= 
tide condition that R«L (where L is the smallest scale length of variation of 
£ and £). But, this does not mean that (6) is valid for any ratio R/L, as is 
discussed in Section 3. (m) The result (6) is usually interpreted to hold 

only for gyrotropic distributions, in which Fjj =Fi _i=0 in (4), (e.g. Rossi anu 
Olbert, 1970, Lee and Fisk, 1981). This too strict condition on F stems from 
imprudent mixing of the <£. v> tensor with the pressure tensor P=N(<p v>-<p><v>). 

( i v ) If F is of the simple form F«(1+cos©) v , the underlined term in T6) exceeds 
the first one in r^j only if Y>.3.56. This represents a much higher anisotropy 
than usually founain Astrophysical applications, because then fully 95% of the 
particles have pitch angles in the forward range O>0>n/2. (v) For first order 

anisotropies, F£j(£>2)=0, the exact guiding centre velocity is form (5): 


or 


< ~9 > 


<v > 

—m 


+ £ x £ + (P/qN) £ x £ 

+ E x b + (pv/3q) V x b 


(7a) 

(7b) 


<V >n 
— it ft 

where P = (N/3)<pv> is the scalar pressure. Because of the error mentioned in 
(m), this is usually quoted as an approximate result, strictly valid for iso= 


tropic distributions only. 


V x b 


i ri 

B _ 


x IB 


(vi 

1 x f 


Since 


(8) 


the last term in (7) is the sum of the usual gradient and curvature drift. 


2. Average Particle Velocity <v> and Nomenclature . From the first moment of 
the Vlasov Equation 3F/3t+v*VF+j)/VpF=0, it is readily shown that in the steady 
state the average particle velocity for a distribution of the form (4) is 

<y.>=<y. || > +Ixb+i|^ : |xV(N<Jp i v i >)+(<p i( V | > -i<P i V i >)l(Vx|) i J (9) 

The underlined term is identical to the one in (6). For first order amsotro= 
pies this reduces to the well-known form 

<v> = <v. > + E x b + 1 / (qN ) b x VP, (10a) 

or in the differential case 

<v> fi = <v. |( > fi + CE x Jb + pv/(3qU) b. x VU (10b) 

where U =p 2 (1/4T)/Fdft is the differential number density in (p,p+dp), and 
C = 1-(1/3U)3/3p(pU) the Compton-Getting factor. The results (10) are standard 
in cosmic ray literature. 

From the definition Xg=x-R it follows that <v.g> = <v> - <d£/dt>. The average 
circulation velocity <dR/dt> is usually called the diamagnetic drift, and to= 
gether with <v>, these two velocities are frequently referred to as eolleetzve 
plasma drzfts'. This is in opposition to the guiding centre drift which is cal = 
led a "single particle drift". The average guiding centre drift seems equally 
collective to us. We also argue that the guiding centre drift is the only drift 
because it is designed to tell you where the distribution goes after the trouble 
some, typically dominant circular motion of each particle has been subtracted. 

In BMW these arguments are extended to the non-steady state (3/3t=0). After 
multiplication of all velocities with qN it is easily shown that for first or= 
der anisotropies 
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i.“ig + i.b + ip“-^g + - x — +9 -^ 9t » • (11) 

where M=-N<Jp^v >b is the diamagnetic moment per unit volume and n=bxN<£> the po= 
larisation vector of the plasma. The total plasma current density is therefore 
the sum of the guiding centre (or free) current density, the diamagnetic (or 
bound) current density and the polarisation current density. When standard 
"single particle" drift analyses are used, the equivalent result for (11) can on= 
ly be written down readily for the perpendicular components (e.g. Parker, 1957). 
Northrop (1961, 1963) also included the parallel components but he needed the in= 
troduction of a distribution function Fg for guiding centres in addition to F, 
and elaborate algebra to do it. By contrast our result is essentially true per 
definition. 

3. Neutral Sheet Drift . Standard drift theory breaks down at neutral sheets 
( B =0 ) because the condition R/L«1 is violated. Our derivations are independent 
of this ratio but they still break down. The reason follows from Figure 1 which 
shows a field B.=B(x)e. z , with B(x) homogeneous in x, and dB/dx any arbitrary large 
positive value. The drift of particles such as (a) and (b) is perfectly well 
described by the last term of (7), independent of the magnitude of dB/dx. How= 
ever (8) shows that at x=0, £xb.-*=°, while Figure 1 suggests that individuals drift 
with finite velocity. Consequently <Vg> cannot be infinite (as it was used by 
Jokipii and Thomas (1981, and references therein)). The obvious reason for this 
breakdown is that the drift velocity at x=0 depends on the phase angle at the 
point of crossing. 

The drift in and around a neutral sheet is readily calculated if B_ is given by 
B(x)=[^2H(x)-1]B e , where B 0 is a constant and H the unit step function at x=0. 
Figure 2 shows the" orbit of a particle with velocity v=(v,0,<t>), projected onto 
the plane 0=tt/2. It last crossed the sheet at point a under phase angle (f>c» and 
will again do so at point b after having completed a projected arc length 
s-2Ro(Tr-cj> c )sin0, where R 0 =p/(qB 0 ). At that time the guiding centre abruptly 
j'umps a distance t =2RoSin0sin<j) C in the direction of -ey. Therefore the component 
of the (time averaged) guiding centre velocity in the direction -ey may be taken 
as 

v gy = f /At = U/s) v sin 0 = v/ (H-<f> c ) ]sm 0 sin 4> c (12) 

From the relation cos<j> c =cos<t>+x/(Rosin0) it follows that of all the particles mo= 
mentarily inside d 3 x, at a distance x from the sheet, only those with 

0i<0<n-0i and 4>i<<()<2ll-<f)x , where sin0i=x/(2R o ) and coscpi=1 -x/ (R o sin0) (13) 

can cross the sheet before completing a full orbit. Therefore the directionally 
averaged drift velocity over an isotropic particle distribution inside d 3 x is 

R-0i 2n-<}> 1 

<v q > = -e y (v/4II) / d0sin 2 0 f dtsin<f> c (0,<j>)/rn-<J> c (0,<}>) 1 (14) 

' “ Si <h 

The x-component of the drift velocity oscillates and does not contribute to the 
average, while the z-component of the average is zero for an isotropic distribu= 
tion. The integral (14) can be evaluated as a series at x=0 to give <v a >=-0,463ey 
while its numerically calculated value in the range -2 Rq<j<<.2R 0 is shown^in Figured. 

This interpretation of guiding centre drift at a neutral sheet differs considera= 
bly from previous ones. Firstly, there is no sheet current density of the form 
J.(x)=J 0 <5(x)ey as is sometimes suggested. In fact the assumed isotropy implies 
<3=0 everywhere, even in a volume element d 3 x at x=0. Secondly, the distribution 
in the entire range -2R 0 <x<2R 0 progresses (drifts) in the direction -e.y, instead 
of the infinite sheet drift current ilq=Jq 0 6(x)ey derived by Jokipn et at. (1977) 
from the V x B/B 2 term in (F). 

As an example we consider the motion of cosmic rays in the interplanetary magnetic 
field. 
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B = B 0 (r 0 /r) 2 £e r - (fir sin 0/V □ P 1 “ 2H(0— rr/2 ) 2 , (15) 

where a is the angular velocity of the Sun, V. is the radial solar wind velocity 
and the function (1-2H) models a flat neutral plane at 0=ir/2. The neutral sheet 
drift (14) will be experienced by cosmic ray particles within an angle A0=2R o /r 
from the plane. If V=400 kms -1 , fi=(2ir/27) day- 1 , r=1 A.U., and B 0 =5//2 nT 
(| B |= 5 nT at Earth), it is readily shown that 

A0 ~ 0.72 (P/P 0 )r(r 2 + r 0 2 )"^ degrees, (16) 

where P is particle rigidity with P 0 =1 GV. Typical cosmic ray protons with kine= 
tic energy T=250 MeV, P=0.75 GV are therefore convected along the neutral plane 
in a region with thickness 2A0~1 degree. This neutral plane drift effect is not 
included in existing drift models of cosmic ray modulation. Jokipn and Thomas 
(1981, and references therein) used the 6-function type drift velocity at the 
neutral plane, derived from the VxB/B 2 term. Potgieter and Moraal (1985) used a 
finite neutral plane drift velocity pattern, similar to that of Figure 2, but not 
based on the fundamental arguments used here. In their model the angular region 
A0 in (16) could l.a. be varied in an empirical and unjustified way, and it is 
significant that cosmic ray observations were better exDlained with finite neu= 
tral plane drift in a region A0£l°, than with a 6-function type drift at 



Figure 2 
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